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ABSTRACT 

We present the results of concurrent X-ray and optical monitoring of the Seyfert 1 
^ ' galaxy Mrk 79 over a period of more than five years. We find that on short to medium 

I time-scales (days to a few tens of days) the 2-10 keV X-ray and optical u and V band 

fluxes are significantly correlated, with a delay between the bands consistent with zero 
OO ' days. We show that most of these variations may be well reproduced by a model where 

I the short-term optical variations originate from reprocessing of X-rays by an optically 

^SJ ■ thick accretion disc. The optical light curves, however, also display long time-scale 

variations over thousands of days, which are not present in the X-ray light curve. 
OO ', These optical variations must originate from an independent variability mechanism 

' and we show that they can be produced by variations in the (geometrically) thin 

disc accretion rate as well as by varying reprocessed fractions through changes in the 
location of the X-ray corona. 
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1 INTRODUCTION 

Active Galactic Nuclei (AGN) emit across almost the en- 
tire electromagnetic spectrum. In most AGN, the bulk of 
the emission is radiated in the UV/optical region in what 
is known as the 'big blue bump'. The blue bump is gen- 
erally accepted to be thermal emission, originating in an 
optically thick accretion disc surrounding the central black 
hole. For a review, see Koratkar & Blaes (1999). AGN have 
been known to exhi bit optical variability s ince shortly after 
their discovery (e.g. lSmith fc Hoffleit|[l96^ ) but the origin of 
the variability is still unclear. The two main possibilities are 
that either the optical variations arise from intrinsic vari- 
ations of the disc emission, perhaps due to accretion rate 
variations, or that the variations arise from reprocessing of 
the variable X-ray emission, boosting the intrinsic thermal 
emission from the disc. 

* E-mail: ebreedt(9astro. soton.ac.uk 



If the optical variations arise from reprocessing of X- 
rays, we expect them to lag behind the X-ray variations. As 
the X-rays originate from close to the central black hole and 
the longer wavelength emission comes from further out in 
the disc, light travel time arguments for a standard optically 
thick disc pr edict that the lag w ill increase with wavelength 
as r cx X^^^ (jCoUier et al.lll999i ). For a typical AGN we ex- 
pect short delays, of order a day or less, between adjacent 
photometric bands. Such wavelength-dependent dela ys have 
been seen in many AG N (Wanders ot al. 1997; Collie r et al. 
1998; 'Kriss et al. 2000 ; Collier ct al. 2001; Oknvanski i et al. 
2003. : Sergeev et al., .2005, ; Doroshenko et al.. ,2006). For a 
detailed discu s sion of the di s c rep rocessing model, see 
ICackett et al] l|2007l '). iGaskelll |2003) suggested that the 
wavelength-dependent delays in the UV-optical may also be 
due to a contribution from reprocessing material much fur- 
ther out in the torus. 

On the other hand, if the optical variations arise from 
variations in the intrinsic disc emission as accretion rate per- 
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turbations propagate inwards, we expect the optical varia- 
tions to precede the X-ray variations. By the same argument, 
the longer wavelength optical emission will also lead those at 
shorter wavelengths. The variations will propagate through 
the disc on a viscous, rather than light-travel time-scale, 
so we expect much longer delays (~months-years) between 
wavebands than in the reprocessing scenario. 

Coordinated X-ray and optical monitoring observations 
provide a method for determining the sign of the delay 
between the optical and X-ray variations and hence for 
determining the dominant process underlying optical vari- 
ability. However, previous studies have yielded varied and 
confusing results, with some sources showing well corre- 
late d emission and othe rs no correlation at all. For exam- 
ple, iMason et aD lj2002h report a 0.14 day lag of the UV 
emission relative to the X-rays in NGC 4051, but on longer 
time- scales [Sh crnmor ct jil. (200^) find a 2.4 day optical 
lead. IPeterson et al.l l|2000i l find a correlation between the 
X-rays and optical on time-scales of a few weeks in the 
same sourc e, but little a g reeme nt on time-scales of days. For 
NGC 3516, iMaoz et all (|20o3) reported a possible 100 day 
X-ray to optical lag after their 1.5 year monitoring of this 
source. However, with the addition of 3.5 years more data 
l|Maoz et al.ll2002 ). they were unable to confirm this lag, or 
measure any other statistically significant correlation at any 
lag in the light curves. For 3C 390.3, the relationship be- 
tween passbands is also complicated: X-ray fiares sometim e 
lead and sometimes follow optical/UV events l|Gaskellll2008h . 
lUttlev et al.l (|2003l ) found the X-ray and optical continuum 
in NGC 5548 to be well correlated on month time-scales, and 
also that the amplitude of the optical variations was larger 
than that o f the X-rays on t h ese lo ng time-scales. On short 
time-scales ISuganuma et al.] (|2006l ) report a clear 1.6 day 
optical lag to the X-ray variations and an optical amplitude 
much smaller than the X-ray amplitude, suggesting that, 
unlike on the longer time-scales, reprocessing is the main 
driver of the optical variability on time-scales of days. In- 
tensive monitoring bv lGaskelll|2006l ). however, fails to show 
this short time-scale correlation. 

Besides simply measuring the optical-X-ray lag, coor- 
dinated X-ray and optical monitoring observations enable 
us to put other useful constraints on the X-ray and opti- 
cal emission processes and on the geometry of the emitting 
regions. For example, the greater fractional variability of 
the X-rays on short (~day) time-scales indicates that, in 
a Compton scattering scenario for the X-ray production, 
at least part of the X-ray variability is not produced by 
variable up-scattered seed photons, but by variations in the 
scattering corona itself. The relative amplitude of optical 
and X-ray variability is also important because if the op- 
tical variations, which represent variations in the bolomet- 
rically dominant blue bump, exceed those in the X-rays, it 
is hard to explain the optical vari ations as purely d ue to 
reprocessing, e.g. as in NGC 554 8 jUttlev et al.l 120031 ) and 
MR2251-178 (|Arevalo et al.ll2008l 'l. In these AGN the large 
amplitude optical variations occur on time-scales of years 
and may represent accretion rate fiuctuations propagating 
through the emission region. 

It has been suggested (e.g. lUttlev et alll2003r ) that the 
strength of the correlation may be dependent on the mass of 
the central black hole. Scaled in terms of the gravitational 
radius, AGN of higher black hole mass or lower accretion 



rate, with cooler discs, will have the thermal optical emitting 
region closer to the black hole than the less massive/high ac- 
cretion rate objects. We may expect stronger correlations in 
the high mass systems, as the accretion rate variations will 
be better connected between the X-ray and optical emitting 
regions. Conversely, in the low mass/high accretion rate sys- 
tems, the very different local viscous time-scales will allow 
the emission regions to vary independently, causing poorer 
correlations on long time-scales. 

In order to try and understand the optical variability 
of AGN, and the relative importance of the two proposed 
variability mechanisms, we have set up a long term mon- 
itoring programme of a sample of AGN of differing mass 
and fractional accr etion rate rh = M/Mgjjj. Previously 
(|Uttlev et al.l l2003h we have reported on observations of 
NGC 5548 \M ~ 7 x l O'^Afa; rh ~ 0.08 (|Peterson et all 
l2004l : IWoo fc Urrvll2002l)] and MR2251-178 [M ~ 10® M©; 
rh ~ 0.04] jArevalo et al.ll2008l ) and here we report on the 
observations of Mrk 79. 

Mrk 79 (UGC 3973, MCG-h08-14-033) is a nearby 
{z = 0.022) Seyfert 1.2 ga laxy with a black hole mass of 
(5.24 ± 1.44) X 10^ Mq (Peterson et al.ll2004h a.nd an aver- 
age M of ~ 6% of its Eddington accretion rate (|Kaspi et al.l 
l2000h . A V band image of the galaxy is shown in Figure [1] 
It is a well-known X-ray source and its o ptical variability 
isweU estabfished (|Peterson et all Il998l : IWebb fc MalkanI 
I2OOOI ). We present the combined optical light curves from 
four independent monitoring programs, together spanning 
more than five years, as well as the X-ray monitoring light 
curve as obtained by the RXTE satellite for the same time. 
Data collection, construction of the light curves and basic 
variability properties are discussed in Section 2. We present 
the cross-correlation analysis in Section 3 and in Section 4 
we show that the medium term fiuctuations are well repro- 
duced by a reprocessing model and that the long time-scale 
trends can be accounted for by allowing the disc and/or 
corona geometry to vary over these time-scales. We conclude 
in Section 5. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 X-ray observations 

We used the Proportional Counter Array (PGA) on board 
RXTE to obtain a long term X-ray light curve of Mrk 79. 
From 2003 March 23 to 2008 January 30, we took snap- 
shots of ^ 1 ks exposure time each, at two day intervals. 
During a period of intensive monitoring, 2005 November 17 
to 2006 January 19, four observations per day were made. 
Only data from PCU 2 were used, as PCUs 1, 3, and 4 were 
regularly switched off during our observations and PCU 
has known problems with high background since the loss of 
its Xenon layer. The data were reduced and analysed with 
FTOOLS v. 6. 4, using standard extraction methods and ac- 
ceptance criteria. The background was calculated from the 
most recent background models which corrects for the recent 
problems with the RXTE SAA history file. The final 2-10 
keV fiuxes were calculated by fitting a power law to the ob- 
served spectra. This allows one to take account of changes in 
the instrumental gain over the duration of the monitoring. 
The X-ray light curve shown in Figure [2] is smoothed 
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Figure 2. 2 — 10 keV X-ray (top panel), u (middle panel) and V band (bottom panel) light curves of Mrk 79. The X-ray light curve 
is smoothed with a 4 point (~8 day) boxcar, to highlight the variations on time-scales similar to that of the optical light curves. The 
X-ray light curve is in units of 10~^^ ergs cm~^ , the optical light curves are in mjy. The u band flux is scaled to a comparison star 
in the field but no colour correction has been applied. V flux calculated as described in the text. No host galaxy contribution has been 
subtracted. 



from six different telescopes here to construct an optical 
light curve spanning more than 5 years. A summary of these 
observations may be found in Table [T] and each is discussed 
briefly below. 




N, 

Figure 1. A combined V band image of the Seyfert galaxy 
Markarian 79 (marked with a red circle) as obtained with the Liv- 
erpool Telescope. The bar at the bottom of the image respesents 
a scale of 1 arcminute. The comparison star used for the optical 
photometry is marked with a blue square, and the summed flux 
from the other two bright stars in the image was used to confirm 
that the comparison star did not vary. 

with a 4 point (approximately 8 day) running average to 
show the variations on time-scales of days more clearly. How- 
ever, the unsmoothed hght curve was used in all calculations. 

2.2 Optical observations 

Since 2002, Mrk 79 has been monitored independently by 
several ground-based observatories. We combine the data 



2.2.1 Liverpool Telescope 

Mrk 79 is part of our ongoing optical monitoring programme 
on the Liverpool Telescope (LT), a 2 metre fully robotic tele- 
scope at the Observatorio del Roque de Los Muchachos of 
the Instituto de A strofisica de Canarias, La Palma, Spain 
l|Steele et al.ll2004[) . The telescope is equipped with an opti- 
cal CCD camera, RATCam, with a 2048x2048 pixel back- 
illuminated detector and a scale of 0.135 arcsec/pixel. Ob- 
servations in the SDSS-u and Bessell-V^ filter bands started 
on 2004 December 2 with a four day sampling period. This 
was increeised in March 2007 to daily monitoring in V and 
every two days in u. The latest measurement in the hght 
curve included here was taken on 2007 May 30. 

The data are overscan-trimmed and flat-fielded by the 
LT reduction pipeline. Bias correction is based on the un- 
derscan region only, as the bias frames do not show any 
significant repeatable structure. Images are binned 2x2, re- 
sulting in a 1024 X 1024 pixel image of the 4.6 arcminute field 
of view. Standard rejection criteria are applied: data from 
nights with poor seeing (> 3") or a very bright moon were 
discarded. The median seeing over the period of monitor- 
ing was 1.3". A few hot pixels and one bad column on the 
CCD were removed by linearly interpolating between adja- 
cent pixels. 

We used the IRAfQ package DAOphot to perform aper- 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
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Table 1. Summary of V band observations 



Observation dates 




Aperture 


No. of V band 


(MJD) 


Instrument and Telescope 


diameter 


observations 


52283.5 - 53987.6 


AP7p camera, 0.7-m Crimean Observatory Telescope 


15" 


142 


52890.6 - 53876.3 


Multicolor Imaging Photometer, MAGNUM Telescope 


8".3 


47 


53100.6 - 54198.8 


0.4-m telescope, Lincoln, Nebraska 


16" 


79 


53341.0 - 54250.9 


RATCam, Liverpool Telescope 


12" 


132 


54063.4 - 54129.3 


0.6-m telescope, Mt. Maidanak Observatory 


16" 


22 


54136.3 - 54237.3 


HawkCam, Faulkes Telescope North 


12" 


39 
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Figure 3. V band light curve of Mrk 79 as obtained with the 
Liverpool Telescope. The ratio of the comparison star to two 
check stars in the field is also shown, showing that the decrease in 
brightness of Mrk 79 is real and not due to a variable comparison 
star. 



ture photometry on the images. We summed the flux inside 
a 12" diameter aperture, centred on the nucleus, and com- 
pared it to a comparison star and 2 chec k stars in the field. 
Comp ariso n star fluxes were taken from iDoroshenko et al.l 
l|2005l ') and lMihov fc Slavcheva-Mihoval (|2007^ . The aperture 
size was selected to include the PSF under seeing conditions 
of 3" or less. In order to exclude host galaxy flux from the 
background measurement, the sky background level was not 
measured in an annulus. Instead, we measured and aver- 
aged the background from three 12" diameter apertures, dis- 
placed in various directions around Mrk 79 and away other 
stars in the field. Aligning all of our images to a reference im- 
age before these measurements were taken, ensured that the 
same sky areas were measured in each image. Photometric 
errors are based on photon statistics and the readout noise 
of the detector. Mrk 79 showed a steady decline in bright- 
ness over the period of monitoring. The V band light curve 
from the LT is shown in Figure [31 relative to the comparison 
stars. 



2.2.2 Faulkes Telescope 

We also carried out a short multicolour monitoring program 
in u,B,V ,R and i on the Faulkes Telescope North (Haleakala, 
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Figure 4. uBVRi light curves of Mrk 79 as obtained with the 
Faulkes Telescope, shown with the corresponding section of the 
X-ray light curve. The X-ray light curve is in units of 10~^^ ergs 
cm^^ s~^ , the optical light curves are in mjy. 



Maui, Hawaii) from 2007 February 5 to May 17. The tele- 
scope and camera are identical in specification and design 
to the Liverpool Telescope and data were reduced and anal- 
ysed in the same way. The sampling is roughly every two 
days, with some gaps due to poor weather conditions. 33 
useful epochs in u, 35 in B, 39 in V, 36 in R and 27 in i 
were obtained. The light curves are shown in Figured 



2.2.3 Crimean Observatory Telescope 

We supplement these light curves with V band data from 
the AGN monitoring program on the 0.7-m telescope 
of the Crime a n As trophysical Observatory, described in 
ISergeev et all l|2005h . The B,V,R,Rl,I filter set on the 
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AP7p camera mounted on this telescope are nonstandard 
filters, though the V filter, from which we include data here, 
matches the standard Bessell-V filter closely. The light curve 
is constructed from aperture photometry through a 15" di- 
ameter aperture centred on the galaxy nucleus, relative to 
a comparison star in the field. It covers the period 2002 
January 9 - 2006 September 9 and has a sampling rate of 
roughly every two days. Data after MJD 53058 (2004 Jan- 
uary 18) from this program was previously unpublished. 



were then used to scale that dataset to the 8". 3 diameter 
aperture MAGNUM light curve. 

The total combined V band light curve spans 1967 days 
(approximately 5 years and 5 months) and contains 461 
epochs. It is shown in Figure [2] along with the X-ray and 
u band light curves. Throughout the rest of this paper, we 
will use this combined light curve when we refer to the V 
band light curve. 



2.2.4 MAGNUM Telescope 

We also include unpublished data from the MAGNUM 

i Multicolor A ctive Galactic Nuclei Monitoring) project 
Yoshii et al.l i2003) . The 2 metre telescope is situated at 
the University of Hawaii's Haleakala Observatory, on the is- 
land of Maui, Hawaii, and monitors a large sample of AGN 
in optical and infrared wavebands. The V band light curve 
included here are the observations made during the period 
2003 September 8 - 2006 May 2. The nuclear flux was mea- 
sured within a 8". 3 diameter aperture, relative to a com- 
parison star in the field. Comparison stars were calibrated 
against photometric standard stars and colour corrections 
applied to the fluxes. Full details of the dat a aqui sition and 
reduction are described in ISuganuma et al.l (|2006l ). 



2.2.5 0.6-m telescope, Mt. Maidanak 

Mrk 79 was also observed with the 0.6-m telescope at the 
Mt. Maidanak Observatory, Uzbekistan from 2006 Novem- 
ber 24 to 2007 January 29. Aperture photometry was per- 
formed on these data using a 16" diameter aperture and a 
background annulus bet ween 24 — 40". The meas urements 
were calibrated using the iDoroshenko et al] (|2005l ) compar- 
ison star magnitudes. 



2.2.6 0.4-m telescope, Lincoln, Nebraska 

Finally, we include V band CCD observations obtained with 
the 0.4-m telescope of the University of Nebraska. Details of 
th ese observations and reduction procedures are as given 
m iKlimek et all (|2004l ). The photometric aperture, sky an- 
nulus, and comparison stars were as for the Mt. Maidanak 
observations. 



2.4 Variability properties 

The flux measured within the aperture contains the nuclear 
flux as well as starlight from the part of the host galaxy 
falling inside the aperture. This host galaxy flux is a constant 
contribution to the total measured flux, so although it does 
not affect the results of the cross-correlation analyis (Section 
[3|, it dilutes the variability of the nucleus. We estimate the 
host galaxy contribution as follows. We obtained an archival 
HST image of Mrk 79 from the online MAST archive, taken 
on 2006 November 8 (MJD 54048) with the Advanced Cam- 
era for Surveys through the F550M filter. This is a medium 
width continuum V band filter, centred on 5580A with a 
FWHM of approximately 540A. We compared the fiux mea- 
sured through an 8". 3 aperture (to match the total V band 
light curve) with the flux measured through a 0".25 aperture 
centered on the nucleus. We measure a 43.5% nuclear contri- 
bution to the total optical emission in the 8". 3 aperture. We 
linearly interpolate our V band light curve to estimate the 
flux at the time of the HST observation and take 56.5% of 
this value to be the host galaxy contribution. This equated 
to a 2.54 mJy correction, which was subtracted from the 
whole V band light curve for the purposes of the variability 
calculations below. No colour correction factor has been ap- 
plied to transform this flux to the flux through the broader 
Bessell-1^ fllter. However, assuming the spectrum to be fiat 
over the width of the fllter, the correction to the flux ratio 
as calculated here, will be very small. 

In order to compare the amount of intrinsic variability 
in each light curve we calculate the error-corrected flux vari- 
ance, as a fraction of the mean flux. This quantity is known 
as the fractional variability. 



if) ' 

where the flux variance has the standard definition 



(2) 



2.3 Relative calibration 



(3) 



In order to account for the differences in the calibration and 
photometry aperture with which the individual light curves 
were constructed, we select the MAGNUM flux calibrated 
light curve as the reference dataset and scale all the other 
light curves to match it. We extracted pairs of measure- 
ments, taken within one day of each other, from overlapping 
light curves, and performed a least-squares fit of the equa- 
tion 



Fi =aF2 + b 



(1) 



to the flux pairs. Fi and F2 are the flux measurements from 
the two datasets and a and b the parameters of the fit. The 
best-fit parameters from each pair of overlapping light curves 



and £ is the sum of the squared measurement errors 



JV 



i = l 



(4) 



The sampling statistics and variability properties of each 
light curve is summarised in Table [2] 

The X-ray and total optical light curves are both highly 
variable. The unsmoothed X-ray light curve contains much 
more variability power at high frequencies than the opti- 
cal light curve, but due to the large amplitude of the long 
term optical variations, there is more variability power in the 
optical light curve when corrected for host galaxy fiux. In 
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Table 2. Sampling and variability characteristics 



Light curve 


Total length (d) 


N 


Ai„edia„ (d) 






X-ray 


1775 


975 


1.89 


0.2733 


10.50 


u (total) 


910 


143 


2.47 


0.3107 


3.81 


V (total) 


1967 


382 


2.00 


0.3447" 


4.20" 


Multicolour monitoring: 












u 


103 


47 


1.75 


0.1432 


1.79 


B 


101 


35 


2.00 


0.0968 


1.63 


V 


103 


51 


1.73 


0.0653 


1.30 


R 


100 


36 


2.00 


0.0504 


1.26 


i 


99 


27 


2.93 


0.0544 


1.23 



" Host galaxy flux subtracted. 



the much shorter Faulkes Telescope colour light curves, the 
fractional variability is seen to decrease with wavelength. 
If the variabiliy is intrinsic to the accretion disc, this be- 
haviour is expected, as the longer wavelength emission orig- 
inate from further out in the disc, where the natural time- 
scales are longer. This is also expected from reprocessing 
models: the disc response function broadens for longer wave- 
lengths, s moothing out the va riability at longer wavelenths 
(see, e.g. ICackett et al.|[2007h . We note, however, that no 
correction for host galaxy flux has been made to these light 
curves, so the variability amplitude we measure is reduced 
by this constant offset. Thus the fractional variability quoted 
in Table [5] should be viewed as a lower limit. 



3 CROSS-CORRELATION ANALYSIS 
3.1 X-ray — optical correlation 

In order to confirm and quantify the apparent correlation 
between the X-ray and optical emission, we use the Discrete 
Corre lation Function (DCF) method of lEdelson fc Kroli"3 
(|l988t ). 

We cross-correlate the unsmoothed, unbinned X-ray 
light curve with the observed u and V band light curves 
and bin the resulting discrete correlation function in uni- 
form lag bins. Throughout this paper we will use a positive 
value on the lag axis to indicate a lag of the longer wave- 
length emission behind the shorter wavelength (i.e. a lag of 
the optical variations behind the X-rays, in this case). 

To account for the ~ 100-day long solar conjunction gaps 
in the optical light curves, we divide them into segments de- 
fined by the observation gaps, and calculate the DCF for 
each segment separately. As the stength of the correlation 
may be different on different time scales, we resample the 
X-ray light curve to a uniform observational frequency by 
selecting only one point every two days from the intesively 
sampled period starting just before MJD 53700. This en- 
sures that this segment does not dominate the combined 
correlation function. The final correlation function is then 
the combination of the individual DCFs, weighted by the 
number of points in each lag bin. The DCFs calculated this 
way, using a bin size of 4 days, are shown in Figure O 

The effect of calculating the DCF in segments like this is 
to implicitly remove the low frequency power (i.e. long term 
trends) in the light curve, leading to an improved estimate of 
the short-term lag. The most common way of removing this 
low frequency power is to subtract a low-order polynomial 
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Figure 5. DCFs between the X-ray and u band (top) and X-ray 
and V band (bottom), binned in 4 day lag bins. The horizontal 
dotted, dash-dot and dashed lines are the median, 95% and 99% 
confidence levels, based on 1000 simulations. The vertical dotted 
line indicates the zero lag position. 



function from the light curve prior to calculating the CCF 
(|Welshlll999l ). We prefer the method of segmentation of the 
light curve, however, as it also eliminates the large observa- 
tion gaps from the DCF calculation. We note, however, that 
removing a linear trend or a higher order polynomial from 
the light curve also result in a correlation peak of greater 
than 99% significance, as discussed below. 

The correlation function peaks at DCFmax = 0.496 for 
the X-ray/1/ band correlation and at DCFmax = 0.654 be- 
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tween the X-rays and u band, with no lag between the op- 
tical bands and the X-rays. The DCF peak occurs at 0^4 
days for V and Ol* days for the u band. (Note that, be- 
cause the DCF was used, these errors will always be a mul- 
tiple of the bin size.) The peak of the correlation function 
is thought to provide the most reliable estimate of the lag 
IWelshlFlQQgl ). however, when using the DCF, time resolu- 
tion on time-scales shorter than the bin size is lost, so the 
peak position is affected by the chosen DCF binning fac- 
tor. We therefore deem the measurement of the centroid lag 
more reliable and stable in this case. We calculate the cen- 
troid as the weighted average of all DCF points greater than 
80% of DCFmax. For the X-ray/V-band correlation, the cen- 
troid lag occurs at days and for the X-ray/u-band 
at -1.95^3:^^ days. The errors quoted are the 1 a errors 
calculated using the Flux Randomisatio n/Ran d om S ubset 
Selection (FR/RSS) method of Peterson et all (119981 ) and 
the centroids are the median centroids from 1000 such se- 
lections. 

As a separate check, we used the z-transform e d dis - 
crete correlation function (ZDCF) of lAlexandeil (|l997l ) 
and the interpolation cross-correlation functio n (ICCF; 
ICaskeU fc Sparkelll986l : IWhite fc PetersonI Il994 ) to calcu- 
late the same correlation. The correlation functions resemble 
those calculated using the DCF closely and the lag results 
are consistent between the three methods. The short lags ob- 
tained from both the peak and the centroid measurements 
suggest that the reprocessor must be located near the X-ray 
source in this object, and not as far out as the dust torus. 

To test the sig nificance of the DCF pea ks found, we 
used the method of iTimmer fc Koenid 1199^ to simulate 
X-ray light curves based on an underlying model power spec- 
trum. We use the single-bend power law parameters, mea- 
sured by Summons et al. (in preparation), as appropriate 
for the power spectrum of Mrk 79. To account for long term 
trends in the data, the simulated light curves were gener- 
ated 10 times longer than the real light curve, and then 
sampled in the same way as the observed X-ray light curve. 
Observational noise was added at the level measured in the 
observed light curve. These uncorrelated, simulated X-ray 
light curves were then each cross-correlated with the real, 
observed optical light curve, to estimate the probability of 
obtaining equal or higher correlation peaks by chance, due 
to the red noise nature of the light curves. We record the 
mean, 95% and 99% extremes of IQQQ simulations and plot 
them in dotted, dot-dashed and dashed lines, respectively, 
in Figure [51 Both the u and the V band peaks reach higher 
than 99% of the simulations, showing that the correlations 
are statistically significant at greater than 99% level. 

3.2 Lags between optical bands 

We also present, in Figure |6] the ZDCF between the 
X,u,B,R and i bands, calculated relative to the V band, 
for the data from the Faulkes Telescope. The u and V band 
light curves were supplemented with data from the Liverpool 
Telescope covering the same period. We chose the V band 
as reference light curve as it is the best sampled. The X-ray 
and V band autocorrelation functions (ACF) are also shown. 
Again a positive value on the lag axis represent a lag of the 
variations in the longer wavelength band behind those of 
the shorter wavelength. The interpolation cross-correlation 
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Figure 6. Cross-correlation functions relative to the V band, 
calculated from the 3 months of Faulkes monitoring. Points with 
error bars are calculated from the ZDCF and the solid line from 
the ICCF. A positive value on the lag axis indicate a lag of the 
long wavelength emission behind the short wavelength. The ver- 
tical dotted lines show the position of zero lag. 

funct ion (ICCF; ICaskeU fc Sparkelll986l : IWhite fc PetersonI 
11994 ) is used as a consistency check. It is plotted as a solid 
line in Figure [6] and is in close agreement with the ZDCF 
results. 

If the optical variations are due to reprocessing of X- 
rays by the accretion disc, standard accretion theory for a 
thin viscous disc allows us to calculate the lags expected 
betw een optical bands. In a standard accretion disc model 
(e.g. iFrank et al] 120021 '). the disc can be approximated by 
annuli at different temperatures, each emitting a blackbody 
spectrum. The radial temperature profile of the disc may 
then be described by 

where Rg ~ GM/(? is the gravitational radius of the 
black hole. Using Planck's Law to calculate the peak emis- 
sion, we can use ((5)l to calculate the radius most of this 
emission is coming from. In the reprocessing model, de- 
lays should exist between different wavebands, due to light 
travel time differences between the illuminating source (as- 
sumed to be the X-rays) a nd the region of the disc emitting 
at the given wavelength l|Collier et al.|[l99^ . Hence vari- 
ations observed at longer wavelengths should follow those 
at shorter wavelengths with a delay r = r/c, where r is 
the distance between the locations of peak emission for 
the re levant wavelengths. Using the reverbera tion mapping 
mass (jPeterson et al.ll2004l ) and accretion rate (|Woo fc Urrvl 
|2002| ) of Mrk 79, and assuming the X-ray source to be a 
point source at height hx ~ GRg above the accretion disc 
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and on its axis of symmetry, we expect a V lag of 1.68 days 
behind the X-rays and lags between the optical bands of u 
vs. V = 0.74, B vs. V = 0.43, V vs. R = 0.64 and V vs. 
i — 0.85 days. 

Both figures |4] and [6] show that the variations in the 
optical bands are well correlated. The centroids of the cor- 
relation function, calculated at 80% of the peak, were found 
to be uvs. V = -0.59+1-92, B vs. V = -O.58l2.29, V vs. 
R = 2.50ti g5 and V vs. i = 1.23l:l ll days. AU except the R 
band are consistent with the expected lags for reprocessing 
to within the 1-a FR/RSS errors. The discrepancy in the R 
band result is likely due to the strong, broad Ha emission 
line contained within the R filter passband. The Ha line 
emission originates in the broad line region, which is further 
from the X-ray source than the region of the disc emitting 
primarily at this wavelength. Hence Ha varies with a longer 
lag than the disc R band emission and increases the average 
lag measured in this band. 

From this short segment of the data, we measure the 

V band emission to lead the X-rays by -^5 days (i.e. X vs. 

V = — 5.25l7 gg days), which is not consistent with the re- 
processing picture. We point out, however, that the error 
bars are overlapping with the X vs. V correlation measure- 
ment found in section |3TT] Since the sampling of these short 
light curves are not significantly better than the total V 
band light curve (see Table [2]), we deem the previous result 
of an optical lag of O.OO+4 qq days behind the X-rays more 
robust and reliable, as it was derived from light curves al- 
most 20 times as long. In short sections of the light curve 
the long term variations (which are present in the optical 
light curve but not the X-rays) can influence the lag mea- 
sured, especially if there are not many peaks or troughs in 
that segment to constrain the lag tightly. In the longer light 
curves, these small shifts average out and a more reliable 
estimate of the true lag is obtained. 



4 MODELLING REPROCESSING 

A possible mechanism for the production of optical fluctu- 
ations is thermal reprocessing of variable X-rays shining on 
optically thick material, possibly the accretion disc. We con- 
structed reprocessed optical light curves, using the observed 
X-ray light curve as input, to compare to the observed op- 
tical light curve. We assume that the X-rays are emitted by 
a compact source above the accretion disc and on its axis 
of symmetry while the optical flux arises thermally from 
an optically thick accretion disc. This geometry of the X- 
ray source may be viewed as a first order description of the 
centre of an X-ray corona above the inner part of the ac- 
cretion disc. More realistic geometries, such as a spherical 
corona of radius hx , will have a small effect on the flux cal- 
culated at the inner parts of the disc, but from further out 
in the disc, where most of the emission we are interested 
in is coming from, the corona will increasingly resemble a 
point source. This simplistic description is thus sufficient 
for our purposes here. We follow ed the prescription detailed 
by JCazanas fc NavakshinI l|200ll 'l to account for light travel 
time effects on the reprocessed light curve, as well as geo- 
metrical considerations on the amount of flux received per 
unit area for different locations in the disc and X-ray source 
heights. We added the impinging X-ray flux to the locally 



dissipated flux expected for a thin disc to calculate the black 
body temperature as a function of time and radius, in order 
to produce flux light curves in different optical bands. 

The model parameters are the height of the X-ray 
source hx, thin disc accretion rate m, inclination angle 
to the observer i and inner truncation radius of the thin 
disc i?in. The mass of the black hole determines the size 
of the gravitational radius and therefore the conversion of 
light crossing time to days. We left the m ass fixed at the 
reverberation-mapped value of 5.24x lO^M© (|Peterson et al.l 
l2004h . We varied Rin and hx between 6 and 50 Rg, where 
the gravitional radius Rg — GM/c^, i between and tv/2 
and m between 0.001 and 0.9 times the Eddington-limit ac- 
cretion rate. A distance of 94.4 Mpc was used to convert 
between luminosities and observed fluxes. No correction for 
reddening or dust extinction has been made. The calculated 
reprocessed light curves were resampled to match the obser- 
vation epochs of the corresponding optical light curves. We 
subtracted the 2.54 mjy galaxy contribution (section I2.4|l 
to the V band flux before fitting the reprocessed light curve 
to the data and subsequently added the galaxy flux to the 
reprocessed light curve for plott ing purposes. 

For a standard thin disc (IShakura fc SunvaevI 119731 ) . 
around a black hole of mass 5.24 x lO^M© and with an 
average accretion rate of ~ 5% of the Eddington limit, 50% 
(100%) of the V band emission is expected to arise from 
within 70 (600) Rg. As the light crossing time for one Rg is 
Rg/c ~250s, an edge-on disc would smooth the reprocessed 
V light curve on time-scales of ~ 2 x 70Rg/c = 0.4 days. 
We note that the inclination angle to the observer is only 
used to calculate the correct light travel times from differ- 
ent locations in the disc to the observer. As light crossing 
time of the optically emitting regions is shorter than the 
usual 2-day sampling rate of the X-rays, the smoothing of 
the fluctuations produced by differential light crossing time 
is small and the reprocessed light curve does not depend 
sensitively on i. 

Figure[7]shows the best-fltting reprocessed V light curve 
together with the observed V band data. The correspond- 
ing parameters are hx = 21.5 i = 24°, m = 0.06 and 
Rin = 37.6Rg . The mean predicted optical flux is very simi- 
lar to the observed flux, the ratio of mean fluxes for the light 
curves shown is 0.993. The normalised transfer functions for 
the V and u bands, calculated using these best-fit param- 
eters, are shown in Figure [S] It assumes an X-ray flare of 
observed flux 4xl0~^^ ergs cm^^ s~^ , illuminating the disc 
for 0.01 days. The flrst peak in the transfer functions corre- 
sponds approximately to the time lag of light reprocessed at 
the inner edge of the thin disc on the side nearest to the ob- 
server, and the secondary peak to the light reflected at Rin 
on the far side. The two peaks merge together for smaller 
inner truncation radii. Varying disc parameters has a small 
effect on the shape of the transfer functions and we note 
that in all cases their median lag is less than two days, i.e. 
smaller than the X-ray sampling. Therefore, the main effect 
of changing model parameters on the shape of the result- 
ing optical light curve is to change the relative amount of 
intrinsic to reprocessed optical flux. 

As can be seen from Figure [T] the reprocessed light 
curve reproduces the medium term fluctuations on the or- 
der of a hundred days moderately well, showing that it is 
energetically feasible to produce these optical fluctuations 
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Figure 7. The fluctuations in the observed V band flux (black 
markers with error bars) cannot be reproduced by reprocessing of 
observed X-rays in a simple scenario where the geometry and disc 
accretion rate remain constant. The solid line shows an X-ray to 
V band reprocessed light curve with disc parameters that best fit 
the observed V band light curve. 
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Figure 8. Disc reprocessing transfer functions for the V band 
(solid line) and u band (dashed line) calculated using the best-fit 
parameters of Figure [7] The first peak corresponds to the time 
lag of light reprocessed at the inner edge of the thin disc on the 
side nearest to the observer, and the second peak to the lag of 
the light reprocessed at the inner disc on the far side. 



through reprocessing. The long term trend in the optical, 
however, is not present in the X-ray band and is not re- 
produced by the reprocessed light curve. This trend in the 
optical suggests that either the geometry changes, affecting 
the amount of reprocessed X-rays or that the thin disc in- 
trinsic luminosity changes, over time-scales of years. To test 
these possibilities, we fitted each segment separately, allow- 
ing only one parameter at the time {Rin, rn or hx) to vary 
between the segments. This setup reproduces roughly the 
behaviour of a change in disc/corona geometry or accretion 
rate on time-scales of years, while all other parameters re- 
main the same. These scenarios, where the long term trends 
are produced by varying one model parameter, will produce 
different u band long-term variations. We therefore include 
the u band light curve in the fit in an attempt to repro- 
duce simultaneously both optical bands with the same set 



of model parameters over the three year segments where all 
three light curves overlap. 

We start by allowing only the disc accretion rate to 
vary from one light curve segment to the next, while keeping 
the other two parameters constant between then segments. 
The best-fitting values of the accretion rate were 11% of 
the Eddington limit for the first segment shown in Figure |9l 
4.8% for the second segment and 0.1% for the last segment, 
while keeping i?in = QRg and hx ~ 6Rg- The resulting re- 
processed u and V light curves are shown in the top panel 
of Figure E] The large range in m is a consequence of the 
small truncation radius used, and represents the case where 
almost all of the optical emission in the lowest flux orbit 
comes from reprocessing. Other solutions, where m does not 
vanish, are possible, but require a larger Rin and results in 
greater contribution from the intrinsic disc emission. The 
middle panel on this figure shows the result of fitting the 
different segments with a variable Rin, keeping the accre- 
tion rate and source height fixed. In this case, the accre- 
tion rate was rh— 7.7% and the source height hx = 7.8Rg. 
The best-fitting truncation radii were Rin = 6Rg (minimum 
value allowed) for the first segment 
second and Rin 
panel in Figure |9] shows the result of a varying source height 
hx- The best-fitting values are Rin = 6Rg, rh — 0.01 and 



Rin = l9.6Rg for the 



hx = 21.37?g for the first segment, hx 



ll.SRg for the 



second and hx ~ QRg for the third. Evidently, varying any 
of these model parameters can reproduce broadly the long 
term trends in the V and u fiux. The details of the repro- 
cessed light curves, however, vary. In particular, varying the 
stable disc accretion rate adds different amounts of constant 
flux to the segments, while the amount of reprocessed flux 
remains approximately the same. Therefore, the short term 
fluctuations, arising from reprocessing, are the same size at 
high and low average fiux levels. The opposite is true when 
varying geometry parameters (source height or truncation 
radius) . In these cases the size of the short term fluctuations 
scale with the mean flux. The average u and V flux levels 
are best reproduced by changing the accretion rate between 
the segments while the worst fit is obtained by varying _Rin. 
Changing this parameter has a much stronger effect on the 
u band than the V band, so the change in u flux is too large 
while the corresponding change in V is not large enough. 
The best overall fit, to the shape and the average fiux lev- 
els is obtained by varying hx, though the quality of the fits 
does not allow us to strictly rule out any of the scenarios. 

We also fitted the multicolour light curves obtained with 
te Faulkes Telescope with the reprocessing model, to com- 
firm that this process can produce the short-term fluctu- 
ations in the different optical bands. We fixed the model 
parameters to the values obtained from fitting the third seg- 
ment of the u and V light curves with either variable hx or 
variable m, as explained above, i.e. Rin = 6Rg, m =1% and 



hj. 



6Rg. As high resolution images exist only for the V 



band, we cannot estimate the galaxy contribution to the fiux 
in the other optical bands, so we left the galaxy fiux as a 
free parameter. The resulting galaxy fiuxes are u — 0.27, 
B = 0.67, V = 2.54, R = 4.08 and i = 5.83 mjy. The re- 
processed light curves for each band are shown by the solid 
lines in Figure [TOl 
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Figure 9. U (blue) and V (red) band light curves and predicted 
reprocessed X-ray flux (solid lines) obtained by changing only one 
parameter in the model between the different segments. 



5 CONCLUSIONS 

We have presented the long-term optical (it and V band) and 
X-ray light curves of Mrk 79, obtained by combining data 
from several monitoring programs of this source. The u band 
light curve is the result of 2.5 years of monitoring with the 
Liverpool Telescope (LT) and the V band light curve was 
constructed from the combined data from six ground based 
telescopes, including the LT. The X-ray observations, which 
were made every ~ 2 days, are from our continuing monitor- 
ing program with the RXTE satellite. We further included 
the results from a 3 month multiwavelength program on the 
Faulkes Telescope in tt, B, V , R and i. 

We summarise our findings as follows: 

(i) The light curves are variable on time-scales as short 
as days (X-rays) and as long as thousands of days (optical) , 
with highly correlated variability on time-scales of tens of 
days or months. 

(ii) Long time-scale (~years) variations are clearly seen 
in the optical light curve, but are not present in the X-ray 
light curve. 

(iii) Although there is more variability in the X-rays than 
in the optical bands on short (~day) time-scales, the frac- 
tional variability of the V band light curve on long time- 
scales is greater when corrected for host galaxy flux, due to 
the large amplitude of these long-term optical variations. 
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Figure 10. Reprocessed light curves for the multicolour intensive 
sampling data. The model parameters were fixed to the best- 
fitting parameters to the V band and the galaxy contribution in 
the other bands was allowed to vary as a free parameter. 



(iv) The lags measured between the X-rays and long-term 
optical bands, are consistent with zero days for the data set 
as a whole. Formally we find that the V band lags behind the 
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X-rays by O.OOI4 qJ days and the u band lags by —1.95^ 
days. 

(v) While some strong events on a time-scale of a few 
days are consistent with the zero days lag, there are also 
time periods where the optical and the X-rays match up 
poorly and events are only seen in one passband. 

(vi) The variations between the short multicolour light 
curves are also well correlated. The measured centroid lags, 
r, are consistent with the r oc A*'^'^ prediction of X-ray re- 
processing by an optically thick accretion disc. However, as 
the multicolour light curves are quite short, with few promi- 
nent peaks and troughs to constrain the lags, the error bars 
on the lag measurements are large and it is not possible 
to constrain the reprocessing model particularly tightly. We 
also cannot exclude the possibility that the measured delays 
across the optical bands are affected by a contribution from 
light reprocessed in the dusty torus further out. 

(vii) Optical variations on time-scales of ~months are 
well reproduced by a reprocessing model, but longer term 
variations (~years) cannot be accounted for in this way. 

In order to reproduce the long time-scale optical trends 
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mentioned above, we are forced to include a second mech- 
anism for the production of optical variations. This mecha- 
nism could be a variation in the geometry of either the pur- 
ported accretion disc or the irradiating X-ray source. Chang- 
ing the height of the X-ray source changes the fraction of the 
disc illuminated and hence the amount of optical emission 
produced by reprocessing. The intrinsic disc flux produced 
by viscous dissipation, is governed by rh and changes in the 
inner truncation radius _Rin affects the amount of optical 
emission produced. We have considered the effect of chang- 
ing hx, rh and -Rin individually, in steps between the three 
sections of the optical light curves (Figure [9]). We find that 
when we vary hx, we are able to reproduce the large and 
small fiuctuations of the observed light curves slightly better 
than when we vary the disc parameters. It is not possible, 
however, at this level of modelling, to rule out the other 
two possibilities. A crude model like this cannot reproduce 
the observed light curves exactly. In reality the parameters 
will probably all vary together, rather than one at a time 
as we have done here, and also gradually, rather than in a 
step-wise fashion. The description of the X-ray source as a 
point source above the centre of the disc, although providing 
a reasonable approximation to an extended hemispherical 
corona, may also not be correct. However, the requirement 
for an additional source of optical variation, in addition to 
X-ray reprocessing, is robust. 

Despite the good correlation on short time-scales, there 
are some short time-scale variations which appear to occur 
only in one passband. This mismatch needs to be addressed 
by future studies and ma y be a conseq uence of anisotropic 
X-rays emission (see, e.g. lGaskellllioO^ ). 

Within the standa rd accretion disc model 
jShakura fc SunvaevI Il973l ). the high black hole mass 
and low accretion rate of Mrk 79 implies that it has a cool 
disc. In such a disc, the X-ray and optical emitting regions 
are closer together, in terms of gravitational radii, than 
in a system of low black hole mass and/or high accretion 
rate. Hence both reprocessing and intrinsic variations are 
expected to contribute to the observed optical variability. 
Reprocessing benefits from the large solid angle subtended 
by the disc in the inner regions and intrinsic variations 
propagate on observable time-scales, because of the shorter 
viscous time-scale at smaller disc radii. We have previ- 
ously reported this behaviour in the quasar MR2251-178 
lArevalo et al.l |2008| ). This system has a very massive 
central black hole [M ~ 10^ Mq) and hence a cool disc, 
with the X-ray and optical emitting regions even closer 
together than in Mrk 79, in terms of gravitational radii. 
When accretion rate fluctuations propagate through the 
disc, the X-ray and optical emitting regions are modulated 
together, leading to a good correlation on long time-scales. 
Reprocessing of X-rays imprints small rapid fiuctuations 
onto the optical light curve, resulting in a well correlated 
behaviour on short time-scales as well. The strong X-ray- 
optical correlation we have reported here is in agreement 
with this picture. 

For Mrk 79 we expect that 50% of the V band emis- 
sion should originate within 70ii3 of the black hole. Taking 
a typical value for the viscosity param eter, a ~ 0.1, in an 
optic ally thick, standard accretion disc (jShakura fc SunvaevI 

0.01, we find that the vis- 
5 years. In a geometrically 



Il973l 'l. with thickness {h/Rf 



thin disc, this time-scale is even longer. This (approximate) 
time-scale, on which we would expect accretion rate pertur- 
bations to propagate inwards to the X-ray emitting region, is 
longer than the time over which we have observed the steep 
decrease in optical flux and hence consistent with the fact 
that we have not yet seen a decrease in long term average 
X-ray luminosity. However the connection between the disc 
and the corona is not well understood and it is not clear 
how accretion rate variations, which may affect the optical 
disc emission, will be translated to the X-ray emission. If 
this interpretation is correct, the X-ray light curve should 
decline signiflcantly over the next 5 or more years. 

In subsequent papers we will compare the optical and 
X-ray variations in a number of other AGN of differing black 
hole mass and accretion rate, and hence of differing disc 
structure, in order to determine how disc structure affects 
both the long and short term optical variability of AGN. 
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